Number fluctuations of cold spatially split bosonic objects 

Kaspar Sakmannjf], Alexej I. Streltsovj^l, Ofir. E. Alon^j, and Lorenz S. Cederbaum 1 ^] 

1 Theoretische Chemie, Physikalisch-Chemisches Institut, Universitat Heidelberg, 
Im Neuenheimer Feld 229, D-69120 Heidelberg, Germany and 
2 Department of Physics, University of Haifa at Oranim, Tivon 36006, Israel 

(Dated: August 17, 2011) 

Abstract 

We investigate the number fluctuations of spatially split many-boson systems employing a theo- 
rem about the maximally and minimally attainable variances of an observable. The number fluc- 
tuations of many-boson systems are given for different numbers of lattice sites and both mean-field 
and many-body wave functions. It is shown which states maximize the particle number fluctuations, 
both in lattices and double- wells. The fragmentation of the states is discussed, and it is shown that 
the number fluctuations of some fragmented states are identical to those of fully condensed states. 
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I. INTRODUCTION 



Ultracold atoms offer the unique possibility to directly compare theoretical predictions 
about many-body physics with experiments. Many-body effects arise due to the interaction 
between particles and the external trapping potential. Of particular interest is the question 
about the nature of the quantum state present in a given system. For example, in double- 
and multi-well trapping geometries, the ground state is either fragmented or condensed, 
depending on the barrier height and the interaction strength jl-7|. For long-range inter- 
actions even the ground state in a single- well trap can be fragmented 8|. Apart from the 
fragmentation of the ground state, fragmentation of Bose-Einstein condensates (BECs) is 



18] . On the theoretical level, frag- 



also known to develop in nonequilibrium dynamics 
mentation manifests itself in the reduced density matrices of the system. While the reduced 
density matrices of a system of bosons themselves are not experimentally accessible, it is 
possible to draw conclusions about them from the measurement of experimentally accessible 
quantities, such as the particle number fluctuations. 



Here, we would like to investigate the number fluctuations of several many-body states 
and their fragmentation. We focus on cold spatially split bosonic objects. Atom number 
fluctuations of fragmented and condensed systems have been investigated intensively both 



theoretically, see e.g. 



15 



19M22| and experimentally 



18 



23 



-28], to name just a 



few. Here, we would like to elucidate the limits that quantum mechanics puts on number 
fluctuations of bosonic mean-field and many-body wave functions, and concentrate on their 
fragmentation. In particular, we find that some fragmented states cannot be distinguished 
by their number fluctuations from fully condensed superfluid states. 



This paper is organized as follows. In Sec. |TT]we discuss a theorem about the maximum 
variance of an observable. In Sec. IIHI some basic definitions needed for the discussion of 
number fluctuations in multi-well traps are introduced. In Sec. IIVI we discuss the number 
fluctuations of two-mode many-boson states. In Sec. |V]we show how the previously obtained 
results can be generalized to BECs in multi-well traps and optical lattices. We summarize 
the results and conclude in Sec. [VH 
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II. MAXIMUM VARIANCE THEOREM 



In this section we will prove a theorem about the maxima and minima of the variance of 
an observable. A more mathematically oriented proof can be found in Ref. {29 1. 

Theorem 1. Let A denote a hermitian operator and ao < . . . < a at a subset of its spectrum. 
Furthermore, let be a wave function and 

l*> = Ec„|* B > (i) 

n=0 

be an expansion of \^) in A's eigenstates with Y<n=o \Cn\ 2 — 1> where C n \^ n ) denotes the 
contribution of all degenerate eigenstates of A to the eigenvalue a n . Then the variance of A 

AA 2 = (A- (^|i|^)) 2 |^) (2) 

takes on its minimum, AA 2 min = 0, for \C m \ 2 = 1, where m G 0, 1, . . . , N is arbitrary. The 
maximum variance, AA 2 max = | (ajv — clo) 2 is obtained for \Cq\ 2 = \Cn\ 2 = 1/2. 

Proof. The normalization constraint X) n=0 \Cn\ 2 = 1 can be used to eliminate one of the 
N + 1 coefficients C n in Eq. (T5]). Without loss of generality, we choose the coefficient Cn 
and write |Cat| 2 = 1 — J2n=o \Cn\ 2 - The variance A A 2 can then be written as 

AA 2 = (^\A 2 \y) - (^\A\q) 2 

N / N 



N I N \ * 

\C n \ a n ~~ I 5^ \C n \ a n J 

n=0 \n=0 / 

N-l /N-l 

^ \C n \ («n — O-n) ~ I X] \C n \ (dn 
n=0 \n=0 



a N ) 



ee AA 2 , (3) 
where AA 2 is a function of 2N independent parameters Co, . . . , CV-i, Co, . . . , C N _ 1 and 

dAA 

dc* 



Y.n=o \C n \ 2 < 1- For an extremum ^Nr = must hold which yields 



JV-l 

C n (a„ - a N ) = 2C n ^ \C m \ 2 {a m - a N ) (4) 

for n — 0, . . . , N — 1. The conditions = yield the complex conjugates of Eqs. (BJ. 

The set of Eqs. (j3J) have the solution C n = for n = 0, . . . , N — 1. In this case |Cv| = 1 
and AA 2 = 0, which is a minimum because AA 2 > 0. Since the choice to eliminate the 



coefficient Cn was arbitrary, any state with \C n \ = 1 for n E 0, . . . , N minimizes AA 2 with 
the value AA 2 min = 0. This concludes the proof of the first part of Theorem [TJ 

Now suppose that d ^ for at least one i E 0, . . . , N — 1, then it follows from Eqs. (J4j) 
that 

N-l 

(di - a N ) = 2 \C m \ 2 {a m - a N ). (5) 

m=0 

In general Eq. (jSJ) can only be fulfilled if |Cj| 2 = 1/2 and C m = for all m E 1, . . . , iV — 1 with 
m 7^ i. It then follows from the normalization condition that |Cat| 2 = 1/2. Since the choices 
i and TV were arbitrary, the maximum variance is obtained for that pair of coefficients 
with |Cj| 2 = \Cj\ 2 = 1/2 that maximizes AA 2 . The maximum variance of A can therefore 
be written as 

AA2 max = \ (ty - a*? (6) 

Since the eigenvalues a n are ordered increasingly, the maximum value of AA 2 max is obtained 
for the choice i — and j = N, i.e. for | C 1 2 = |Cat| 2 = 1/2. This concludes the proof of 
Theorem [TJ □ 



III. DEFINITIONS 



In this section we briefly recall some definitions that will be useful for what follows. Let 
the operators b\ and % denote the operators that create and annihilate a boson in the orbital 
(pi and fulfill the usual bosonic commutation relations [hi, K] = S^. The number operator of 
the orbital 0, is hi = b\i>i. For a permanent in which iV bosons reside in s orbitals with rij 
bosons in the orbital 0« we use the shorthand notation 

\r\2 n \...,s n °) = ° 1 / , 2 ° s 10). (7) 
Vni!n 2 ! • • -n s \ 

The most general wave function of iV identical bosons expanded in s orbitals then reads 

N 

m= E c ni ... ns \r\2 n \...,s n *), (8) 

ni,...,n s =0 

with rii + . . . + n s = N. Furthermore, we write 



if(x) = Y,Ui(x) (9) 



for the bosonic field operator and for an iV-boson wave function. Then p(x) = W(x)if?(x) 
is the operator of the single-particle density p(x) = (p(x)). The first-order reduced density 
matrix (RDM) is defined as 

p {1 \x\x') = (V\&(a/)i?(x)\%) 

= Y,T>>i 1) <*i 1 \x)<*i 1) *W) (10) 
i 

and has eigenfunctions oq (x) and eigenvalues which are known as natural orbitals and 
natural occupation numbers, respectively. Explicitly, 

' dx' p {1 \x\x')af\x') = n^a^\x) (11) 



holds, where . is assumed and Ei nY> = N. If an eigenvalue n. 



exists the system is said to be condensed 



BEC is said to be fragmented 



(i) 



30 



O(N) 



31] , see also 



3 . If there is more than one such eigenvalue, the 
5K7|. The density fluctuations are given by 

Ap 2 (x) = (p(x) 2 ) - (p(x)) 2 . (12) 

In practice, Ap 2 (x) must be integrated over some finite region of space. If each orbital 4>i(x) 
is localized around x = X{ and has little overlap with other orbitals, the integral of Ap 2 (x) 
over a region of space where 4>i(x) is not negligible can be approximated by 

An? = (nf> - (n>) 2 . (13) 



The quantities An 2 are known as number fluctuations and will be discussed in the following. 



IV. TWO-MODE STATES 

Let us now investigate the number fluctuations and the fragmentation of some particular 
many-boson states constructed either from two localized modes, denoted </>£, and 0^ with 
<Pl(z) — <Pr{— x ) or their gerade and ungerade combinations, denoted <j) g {x) = ^[<I>r{x) + 
4>l(x)] and (j) u {x) = -^[(Pr(x) — <Pl(x)). The particle number operator of the orbital 4>R can 
then be written as 

n R = b ] R b R = ^(n g + fi u + b\b u + P u b g ^ . (14) 
Furthermore, if only two modes are available it follows from hi = N — ur that 

An 2 L = An% (15) 
irrespective of the quantum state. We will therefore drop the site index in this section. 
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A. Many-body states 



In recent theoretical work based on the time-dependent many-body Schrodinger equation 
it was shown that superpositions of macroscopic quantum states can be created by scattering 
an attractively interacting BEC from a barrier 131, |32|. The resulting state is known as a 
colon and has two dominant contributions in the basis of left and right localized orbitals. 
We idealize this caton state here by 

1 



I***) = (\L N ) + \R N }) , (16) 

which is also known as a NOON state, since it can be written as ^= (|iV, 0) + |0, N)) using 
the conventional number state notation. For the state \^ C at) we find that the number 
fluctuations are given by 

An lat = N 2 /4. (17) 

Since |0, N) and \N, 0) are the eigenstates of ur and hi, corresponding to their minimal and 
maximal eigenvalues, it follows from Theorem [1] that the state \^ C at) is a state that maximizes 
the variance of hn and hi. Note that any state of the form -^[\L N ) + exp(i9)\R N )} leads 
to the same number fluctuations An 2 = N 2 /4. The first-order RDM of such states has two 
macroscopic eigenvalues nl = r4 — N/2, and thus the caton is a fragmented BEC that 
maximizes the number fluctuations. Thus, a measurement of the number fluctuations is 
insensitive to the relative phase 9 between |JV, 0) and |0,iV). 

This result should be compared to that of a caton state in the basis of the gerade and 
ungerade orbitals 4> g (x) and u (x) 

= ^ (\g N ) + \u N }) , (is) 

which has exactly the same set of eigenvalues of the first-order RDM Tip = = N/2, but 
much smaller number fluctuations which are given by 

AnJ /uoai = l -N. (19) 

Note that any state of the form [\g N ) + exp(i9)\u N )} has the same fragmentation and 
number fluctuations as \^ g / UC at)- Moreover, it is easy to see that also any state of the form 
cos(9)\g N ) + sm(9)\u N ) has the same number fluctuations, An 2 = N/4. So far, no caton 
states have been reported in experiments. Equations (fT7j) and (Tl9|) and the considerations 
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above clearly show that caton states cannot be characterized uniquely by their number 
fluctuations, or their fragmentation ratios alone. 

B. Mean-field states 

Spatially split mean-field states that have received a lot of attention are the soliton train 

states 

l*i) = l^>, \Kt) = \u N ) (20) 

that describe spatially split, fully condensed BECs, i.e. condensates with = N. Soliton 
trains appear in the context of attractively interacting BECs within the framework of Gross- 
Pitaevskii theory. Interestingly, one finds for their number fluctuations 

A< = \n (21) 

which is exactly the same result as for the \^g/ uca t) state. Similar to the case of caton 
states discussed above, also the state \u N ) leads to An 2 = jN. Thus, a measurement of the 
number fluctuations alone does not allow to distinguish between the states \^ g / uca t), \^ft) 
and \^st)- A simultaneous measurement of number fluctuations and fragmentation would 
be necessary to narrow down the number of possible states that the system was in. 

Let us now turn to more general mean-field states. To this end we define parameterized 
two-mode operators 

a x (9) = cos(9)b L + sin(9)b R , 

02(B) = - sin(9)b L + cos(9)b R (22) 

which can annihilate bosons either in localized or delocalized orbitals depending on the value 
of 9, and compute the number fluctuations of the general mean-field state \^mf) that can 
be constructed from the operators a\ and a\: 

\y MF ) = \a1(9),a%- n (9)). (23) 

The number fluctuations of \^mf) are given by 

sin 2 (29). (24) 



An 2 - 

Lxn MF — 



N n(N - n) 

T + 2 
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The maximum of the number fluctuations An 2 MF when considered as a function of 9 and n 



is given by 



.2 N 2 N 
max An MF = — + — (25) 



which is obtained for n = N/2 and 9 = 7r/4. For these values of n and 9 the wave function 
\^mf) becomes 

W = l^.^)- (26) 

The state l^/mg) is a so called fragmenton state [Hj. For attractively interacting BECs it 
was recently shown that soliton train states can quickly loose their coherence and become 
spatially split, fragmented objects, like the fragmenton state 17fl. Just like the two caton 



states discussed above, fragmentons are fragmented BECs. The state \^f rag ) is two- fold 
fragmented with rii = t4 = N/2. Interestingly, the number fluctuations of the fragmenton 
frag) has contributions oc N and oc iV 2 , see Eq. ( 1251) . The minima of An 2 MF are obtained 
for 9 = 0,7r/2,... with n G 0,...,N arbitrary. The corresponding states are known as 
number states or Fock states 

\* nwn ) = \L n ,R N - n ) (27) 

and have zero number fluctuations An 2 = 0. The fragmentation of number states depends 
on the number of particles in each localized orbital, and is given by n x = n, = N — n. 
This concludes our discussion of two-mode systems. 



V. LATTICE STATES 

A. General lattice states 

We will now generalize the discussion to lattices with s sites, denoted i = 1,2, ...,s 
and corresponding localized orbitals 4>i(x). First, we will prove that the maximum num- 
ber fluctuations in an s-site lattice are identical to those of a two-mode system, namely 
max An 2 = N 2 /4. For simplicity, we begin with a lattice of s = 3 sites. The ansatz wave 
function, Eq. (jSJ), then reads 

N 

|* 3 ) = E a™|l n \2 n2 ,3 n3 >, (28) 

ni,n 2 ,n 3 =0 

8 



with m + n 2 + n 3 = N and Eni,n 2 ,n 3 \C ni n 2 n 3 \ 2 = 1- We define 



N—n 



|Cni| — ^ ] | Cni,Ti2,jV— ni— 712 1 (29) 
ri2=0 

and note that Y^ni=o l^nil 2 = !■ The variance An 2 of (ni) can then be written as 

An\ = ((ni - (r^)) 2 ) 

TV / JV \ 2 

ni=0 \m=0 J 

Analogous to our proof of Theorem (TJ it follows that for a three-site lattice the maximum 
number fluctuations are An 2 = iV 2 /4, obtained for |C ni= o| 2 = |C ni= 7v| 2 = 1/2. Similarly, the 
minimum number fluctuations are An 2 = 0, obtained for \C ni | 2 = 1 for any n\ G 1, . . . , N. 
Since the choice of the lattice site i — 1 was arbitrary, the minimum and maximum number 
fluctuations at any lattice site i — 1, 2, 3 are Anf = and An 2 = iV 2 /4, respectively. Thus, 
we recover the same values for the minimum and maximum number fluctuations as in the 
case of two lattice sites, see Sec. IIVI For a lattice with s sites the same reasoning applies if 
|C ni | 2 is redefined as the sum over all absolute value squares of coefficients with exactly n\ 
bosons at lattice site % = 1 [see Eq. (I34p below]. Thus, we find 

iV 2 

max An 2 = — minAn 2 = 0, (31) 

for the absolute maximum and minimum number fluctuations for lattices with s sites. 

In the present calculation no assumption was made about the symmetry of the wave 
function. Hence, states that maximize the number fluctuations An 2 will generally have 
different number fluctuations at different lattice sites. This can easily be seen by noticing 
that, e.g., the state ^(| AT, 0,0) + |0, 0,iV)) has number fluctuations An\ = N 2 /A at the 
first, but An 2 , = at the second lattice site. States that maximize number fluctuations and 
possess the symmetry of the lattice will be treated next. 



B. Symmetry restricted lattice states 

We will now require that all lattice sites be equivalent with mean occupation {hi) = N/s. 
Since all sites are assumed to be equivalent we will drop the site index from now on. As 
before, we begin with s = 3 lattice sites. It is easy to see that the three-site lattice caton 
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state 

' I JV\ , loiV\ , \oN 



1 /, 



\y cat _ 3 ) = — (\l«) + \2 N ) + \3 N )) (32) 
has mean occupation (n) = N/3 and number fluctuations 

AnL-3 = In 2 (33) 

for all three sites. Its number fluctuations are slightly less than the maximal possible value 
An 2 = N 2 /4, and we will now show that these are also the maximum number fluctuations 
under the constraint of equivalent sites. 

As an ansatz for the wave function on the lattice we use Eq. (JSJ. Let us focus again on 
the number fluctuations on one, say the first, of the s equivalent lattice sites and define the 
quantities 

N-ni N—ni—...—n s -i 

\Cni\ = E ■ ' ' E |Cm,...,Ar-ni-...-n s _i| • (34) 

"2=0 n s _i=0 

The requirement of mean occupation N/ s on all lattice sites can be written as 

N N 

E \C n fn t -- = (35) 

n,=0 S 

for i = 1, . . . , s. Using the equivalence of all sites, we can focus on the first lattice site, and 
after dropping the site index the problem reduces to finding the extremum of the functional 



N / N \ 

r[{C* n }, {C n }} = E \C n \ 2 n 2 - E l^n| 2 n 

n=0 \n=0 / 

( N — N \ 

-V E l^nl 71 

\n=0 S / 



N 

■ \ 

(36) 



where the normalization Y,n=o \Cn\ 2 = 1 is used. This normalization constraint can be used 
to eliminate |CV| 2 in Eq. (136]) . giving 



N-l /N-l \ 2 

r = E \C n \ 2 (N - nf - E l^nl 2 ^ - n) ) + 

n=0 Vn=0 / 

/N-l 



A*( E \C n \ 2 (N-n)-—N). (37) 



\n=0 

For an extremum dr/dC* = must hold, i.e. 







(N - n) 2 - 2(N - n) ( - — -N ) + fi(N - n) 



s 



C n (38) 



for n = 0, . . . , N — 1. If C n = for n = 0, . . . , N — 1, Eqs. (138|) are satisfied, but it follows 
from the normalization that |Cat| 2 = |Cjv,o,...,o| 2 = 1, i>e. the ansatz wave function, Eq. (jSJ), 
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reduces to \N, 0, . . . , 0). Not all sites are equivalent in the state \N, 0, . . . , 0) and therefore 
there is no solution with C n = for n = 0, . . . , N — 1. Thus at least one \C n \ 2 must be 
nonzero for n = 1, . . . , N — 1. Assuming one particular nonzero C n , Eq. (138j) puts the 
constraint 

s — 2 

H = + n ( 39 ) 

for each value of n G 1, N — 1 on \i. Obviously, this constraint can only be satisfied 
for at most one n. Thus, solutions to Eqs. f l38l) must be of the form C n ^ and CV ^ 0, 
and the normalization constraint becomes \C n \ 2 + |Cjv| 2 = 1. Likewise, the requirement of 
mean occupation N/s reads — = |C n | 2 ra + |Cjv| 2 A^. The two conditions can be combined to 
express |CV| 2 and \C n \ 2 as 

IN - sn 



C 



N\ 



1 2 



s N — n 
= j _ 1 N-^ 
s N — n 

which in turn can be used to express An 2 = \C n \ 2 n 2 + |Cat| 2 A^ 2 — (|C n | 2 ?7, + \C N \ 2 N) 2 after 



some algebra as 



v 2 s(s — 1)N s — 1 ~ 

An 2 = -n v ^ 2 ; + — —N - (41) 



s 

The maximum of An 2 as a function of n is obtained for n = with 



s — 1 

max An 2 = — ^-N 2 . (42) 



,S' 2 

Substituting n = in Eqs. (jl"0|) we find that the maximum particle number fluctuations for 
states with equivalent lattice sites are obtained for 

l^l 2 = — |CV| 2 = -• (43) 
s s 



Let us now return to the three-site lattice caton state, given in Eq. (1321) . By setting 
s = 3 in Eq. (l4*2|) and comparing the result to Eq. ( 1331) . we find that l^cat-s) is a state 
that maximizes the particle number fluctuations under the constraint that all three sites are 
equivalent. Note that also states with nonzero relative phases between the components of 
l^cats) would lead to the same number fluctuations, but the sites would generally not be 
equivalent then. More generally, we find for the s-site caton state 

1 



\y ca t~s) = ^= s {\l N } + \2 N ) + .-- + \s N }) (44) 
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that the number fluctuations are given by 

An lat-s = ~~2~-^ 2 - (45) 

This means that s-site caton states maximize the number fluctuations under the constraint 
that all sites are equivalent, see Eq. (142 j) . 

Let us now discuss s-site mean-field states. In case that N/s is integer, it is easy to see 
that the Mott insulating state 

\t! MI ) = \l N ' s ,2 N l s ,...,s N l s ) (46) 

is a lattice state with equivalent sites that minimizes the number fluctuations with An 2 MI = 0. 
Therefore, the complete range of number fluctuations under the constraint of equivalent sites 
is 

< An 2 < S -^N 2 . (47) 
Lastly, we define ¥ g = ^(b\ + . . . + t>l) and discuss the superfluid lattice state 

I*./) = \9 N ). (48) 

The ground state of noninteracting bosons in a lattice potential is of this form and its number 
fluctuations are given by 

An 2 f = N^=±. (49) 

By comparison with Eq. (I47j) it becomes clear that the superfluid state \^ s f) is about in the 
middle of the range of possible number fluctuations. It is fully condensed and hence its first 
order RDM has only one macroscopic eigenvalue, = N, i.e. there is no fragmentation. 
The state \^ s f) is by far the most intensively studied state and concludes our investigation 
here. 

VI. CONCLUSIONS 

We have studied the number fluctuations and the fragmentation of various many-boson 
states, focusing on ultracold spatially split systems. Number fluctuations are a key quantity 
in determining the state of a quantum system. We have shown that there is a great inde- 
terminacy if number fluctuations are considered alone. Additional observables will have to 
be considered to allow for conclusive results, e.g. the fragmentation. For an overview of the 
obtained results please see Table HI 
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Object 


# sites 




An 2 


max An 2 


Caton 


2 


N/2 


N 2 /4 




g/u Caton 


2 


N/2 


N/4 








N 2 /4 


Soliton trains 


2 


N 


N/4 


Fragmenton 


2 


N/2 


N 2 /8 + N/4 




Lattice caton 

Mott-insulator 

Superfluid 


s 
s 
s 


N/s 
N/s 
N 


N 2 (s-l)/s 2 


N(s- l)/s 2 


N 2 {s-l)/s 2 



Table I. Number fluctuations and fragmentation of different spatially split bosonic objects. Given 
are the number of sites over which the object is distributed, the largest eigenvalue of the first-order 
reduced density matrix n[\ the number fluctuations An 2 and the maximally obtainable number 
fluctuations max An 2 . Only objects for which all sites are equivalent are shown. 
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